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Abstract 

The C r dependence problem of multiple Dirichlet eigenvalues on domains is discussed for 
elliptic operators by regarding smooth one-parameter families of C 1 perturbations of domains in 
M n . As applications of our main theorem (Theorem 1), we provide a fairly complete description 
for all eigenvalues of the Laplace operator on disks and squares in M 2 and also for its second 
eigenvalue on balls in M" for any n > 3. The central tool used in our proof is a degenerate 
implicit function theorem on Banach spaces of independent interest. 


1 Introduction and statement 

A number of eigenvalue problems associated to the Laplace operator have been widely investi¬ 
gated since the famous Rayleigh’s book [16j. We refer for instance to the Henrot’s [9] book for 
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details on background material. 

The present paper deals with the dependence problem of Dirichlet eigenvalues for elliptic 
operators with respect to perturbations of the domain of embedding type. This is a classical 
problem in the literature which has been addressed at different times, mainly regarding alge¬ 
braically simple eigenvalues, namely eigenvalues with algebraic multiplicity equal to 1. For this 
case we refer among other works to m, 0, 0 , 0 , 1 , 0 , PI. PI and 0 . 

In the case that the eigenvalues are algebraically multiple, much less results on dependence 
are known. Indeed, one does not hope in general that multiple eigenvalues depend differentiably 
on smooth perturbations of the domain. Examples can easily be constructed even in finite 
dimension (see for example page 37 of 0). Nonetheless, weaker differentiability of eigenvalues 
with respect to differentiable deformations of C 2 domains has been investigated in the multiple 
situation by Cox, namely in the sub differentiability (or Lipschitz) sense (see Theorem 1 of [5J), 
and by Rousselet m and Munnier jT33 , in the lateral directional derivative sense (see Theorem 
2.5.8 of 0). On the other hand, the C r dependence of multiple eigenvalues upon smooth one- 
parameter families of perturbations of the domain seems to be unknown even for r = 1 and 
smooth perturbations. Important contributions can also be found in the classical Kato’s book 
m where are provided expansions of eigenvalues considering one-parameter perturbations in 
the analytical context. For an excellent overview on dependence of eigenvalues with respect to 
the domain, among other interesting problems, we also refer to the Henrot’s [S3 book. 

According to the works of Uhlenbeck [181 EH] and Pereira [T3j, it is interesting to note that 
algebraic simplicity of eigenvalues of the Laplace operator is a generic property. Precisely, 
given a bounded open subset 12 of M n and a number £ > 0 , there exists a diffeomorphism 
99 : 12 —» <^(12) C M n such that ||</? — Id\\c3 < £ and 99 (D) has the following property: all 
eigenvalues of the Laplace operator on 99 ( 12 ) are algebraically simple. Nevertheless, eigenvalues 
of symmetric domains (e.g. balls and cubes in M n ), except the first one, have in general algebraic 
multiplicity greater than 1 (see ID- 

More particularly and surprisingly, the following question is open for any r > 1: 

Does the Dirichlet spectrum of the Laplace operator on balls in M n vary C r smoothly upon 
one-parameter perturbations of C 1 class of the domain for any dimension n > 2 ? 

By using an appropriate degenerate implicit function theorem on Banach spaces, we es¬ 
tablish a result on C r dependence of Dirichlet eigenvalues of elliptic operators in the multiple 
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case (Theorem 1). As a byproduct, we give a fairly complete answer, in an almost everywhere 
sense, to the problem stated above in dimension n = 2 on disks and squares and also for the 
second eigenvalue on balls in M n for any dimension n > 3. As we shall see, perturbations of the 
domain (which can be a certain non-smooth domain) mean images by diffeomorphic C 1 maps 
and deformations of these perturbations are considered C r+1 with respect to the parameter. 

Before we go further and state our main theorem, a suitable framework of the problem of 
interest should first be introduced. 

The eigenvalue problem for the Laplace operator under Dirichlet boundary condition is 
given by 

{ —A u = Ao u in fl, 

( 1 ) 

u — 0 on dd, 

where 9 denotes a bounded open subset of M n , not necessarily smooth, n > 2 and Ajj 6 K is a 
eigenvalue of (DO). 

Let h2o be a bounded open subset of M". Roughly, let Q t be a smooth one-parameter 
perturbation of h2 0 around t — 0. We are interested in studying the C r regularity of the 
map t i—^ An t around t — 0 in the case that Aq 0 is an algebraically multiple eigenvalue. The 
techniques often used in the context of algebraically simple eigenvalues base on bifurcation 
theory, implicit function theorem, perturbation theory of operators and transversality theorem. 
Such tools do not apply directly to the multiple case, so instead we approach the problem 
through a degenerate implicit function theorem on Banach spaces to be proved in the Section 
2 (Theorem [2]). The idea is simple and consist in constructing a suitable map with invertible 
derivative by using the first and second derivatives of the map for which we seek an implicit 
function. The method may vaguely resemble a Lyapunov-Schmidt reduction process (to finite 
dimension) but no Fredholm condition is required and also there is no additional nonlinear 
equation in our assumptions. These characteristics give Theorem [2] an independent interest 
and possibly other interesting applications. 

We now make precise the meaning of the expression “smooth one-parameter perturbation of 
Oo” based on the Uhlcnbeck’s idea [T9] of parameterizing domains as images of a fixed domain 
via diffeomorphisms. The advantage is that the collection of such diffeomorphisms is a subset 
of a certain normed vector space and, therefore, has differentiable structure. 

Let fio be a bounded open subset of M n . Denote by X the Banach space of C 1 maps 
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p : h2 0 —> M n , which extend continuously up to the boundary of h2 0 as well as its derivatives of 
first order, endowed with the usual norm 


IMU : = max{||</?(z)||, \\Dp(x)\\} . 

xGQo 

Let i? 1 (h2o) be the collection of maps p G X such that p is a diffeomorphism on its image. 
Note that the inclusion function 1 q 0 belongs to £ ,1 (h2o)- 

For a bounded open subset of we denote by A(—A, h2) C M x the set of couples 

(Aq,w) satisfying (JTj) where u is nonzero. 

Our main theorem states that 

Theorem 1. Let h2o be a bounded open subset of M n satisfying the divergence theorem for an 
unit normal vector field v defined on dLl® in almost everywhere and oriented outwards Oo- Let 
p : M —> £ ,1 (h2 0 ) be a C r+1 curve, r > 1, such that t^(0) = 1 q 0 . Denote Tp = p'(0). Let A 0 be a 
multiple eigenvalue of the Laplace operator on Oo - Assume that its eigenspace is generated by 
a basis {w 0 i Ui, ■ ■ ■, Uk} C Hi (h2 0 ) with k > 1. Consider the (k + 1) x (k + 1) real matrices A 
and B with elements, respectively. 


and 


Aij ( a, . V. j ) <p . 


dui duj _ 

p ■ v db 


Idn 0 dv dv 


L^ij {d,■ uf )],2 . I UiUj dx . 

J fto 

Assume that the polynomial y(s) = det(A — sB ) has a simple real zero p. Assume also that 
the eigenfunctions Ui are of C 1 class up to the boundary. Then, there exist C r functions t 6 
(—£,£) •—y A(f) and t e (—e, e) i —y u(t ) such that A(0) = Ao, w(0) G (uq, u\, ..., Uk) and 
(A (t),u(t)) G A(—A, (/?i(h2 0 )) fort G (—£,£). Moreover, we have A'(0) = —p. 

In the above theorem, the assumption that eigenfunctions corresponding to Ao belong to 
C ,1 (Oo) can be relaxed by assuming a L 2 integrability of trace of their gradient. This is a case 
for a wide class of non-smooth domains. 

The matrix A has already appeared in Theorem 3.2 of PU (see also Theorem 2.5.8 of 0) 
when the basis {u 0 ,ui,... ,Uk} is L 2 -orthonormal, where its eigenvalues are proved to be the 
possible values of the lateral directional derivative of A(f) at t = 0. The statement of Theorem 
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Q] requires the zero /i of y(s) to be simple so that no other branch of eigenvalues A(t) has the 
same derivative at 0. This key assumption allow us to prove the C r smoothness of A(f) around 
t = 0 with the aid of the perturbative tool of the next section. 

The paper is organized into three sections. Section 2 is devoted to the proof of our main 
tool, the degenerate implicit function theorem on Banach spaces. In Section 3 we present an 
extension of Theorem 1 to more general elliptic operators and provide its proof. Finally, in 
Section 4 we solve the proposed problem on disks and squares and also for the second eigenvalue 
on balls in M n for almost every perturbation ip of the domain. 

2 An abstract perturbative theorem 

This section is devoted to the proof of the following degenerate implicit function theorem. For 
the concepts of infinite dimensional submanifold and other definitions from infinite dimensional 
differential geometry, we refer to [T] or m 

Theorem 2. Let X and Y be Banach spaces and F : M x X —y Y be a C r+l map , r > 1, 
and M C X be a C r+l submanifold with complemented tangent space. Let p be a fixed point on 
M and set q = F(0,p) and L = F x (0,p) := D x F(0,p). Assume the following conditions are 
satisfied: 

(a) F( 0,x) = q for all x G M; 

(b) Ker(L) = T P M; 

(c) F t (0,p) := -D ffi F(0,p) = L(v) for some v G X; 

(d) The composition 

T P M Y coker (L) 

is a vector space isomorphism, where G = D(F t — F v )(0,p) and i r is the canonical projec¬ 
tion. 

Then, there exists a C r curve x : (— e, e) — * X defined in a neighborhood of 0 such that x(0) = p 
and F(t,x(t)) = q for all t G (—£,£). Moreover, we have a/(0) + v G T P M. 
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Remark 1. The co-kernel defined as the quotient coker(L) = — needs not in general to be 
a Banach space but just a vector space. In our applications, Im(L) is finite-codimentional and 
thus, complemented. 

Proof of Theorem 2. By direct computation, it is easy to see that all properties are preserved 
by C r+l changes of coordinates. In other words, if 'h : X —> X is a local C r+1 diffeomorphism, 
then the function F(t,x) = /(t, v h _1 (x)) also satisfies the conditions (a)-(d). In fact, if M = 
'I'(M), p = 'k(p), q = q, A = D^(p),L = DF(0,p),v = A(v) and G = D(F t — F c )(0,p), then 
L o A = L and GoA = G + Lo D (\o A, so that the claim follows. 

By taking a suitable coordinate neighborhood at p E U C X adapted to M, we may assume 
that there is a decomposition X = Xi © X 2 such that M HU = ({0} x X 2 ) D U. Moreover, we 
assume that q — 0. By the assumptions (a) and (c), F : R x X x x X 2 —> Y satishes 

F(0, 0, x 2 ) = 0 

for all x 2 G X 2 around 0 and 


*1(0,0,0) = L(v). 

In addition, by (b), the partial derivative of F in the X,-direction at (0,0,0), denoted by 
L Xl : Xi —> Y, is injective since Ker(L) = {0} x X 2 . 

Let w £ Xi be the hrst coordinate of v, so L Xl {w ) = L(v). Now consider the functions 

H(t, xi, x 2 ) = F(t, t(x i - w),x 2 ) 

and 


J(t,xi,x 2 ) = -H(t,xi,x 2 ) = / Hfist, xi, x 2 ) ds 

I Jo 

which is an integral of C r functions. By Example 2.4.16 (Differentiating Under the Integral) in 
[T], the function J is C r . 

On the other hand, 


J(0,xi,x 2 ) = H t { 0,xi,x 2 ) = DF( 0,0,x 2 )[l,xi - uyO] 


and 
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DJ{ 0, 0, 0)[0, x u x 2 } = DF{ 0, 0, 0)[0, x u 0] + D 2 F( 0, 0, 0)[(1, -w, 0), (0, 0, x 2 )\ 


= : L Al [x i] + G[x 2 ] ■ 

Note that the bounded operator L Xl [xi]+G[x 2 ] is a vector space isomorphism by the assump¬ 
tion (d), thus a Banach space isomorphism. In addition, J(0,0,0) = DF(0, 0, 0)[1, — w, 0] = 
F t ( 0, 0, 0) — L Xl [w] = 0. So, the implicit function theorem applied to J provides e > 0 and C r 
functions X\(t) and x 2 (t) such that 

0 = J(t,£i(t),x 2 (t)) = jF(t,t(xi(t) - w),x 2 (t)) 

for all t G (— e, e) with t ^ 0. Therefore, the result follows by taking x(t) = (xi(t), x 2 (t)), where 
xi(t) = t(xi(t) — w) and x 2 (t) = x 2 {t). Indeed, the above equality yields F(t,x(t)) = 0 for all 
t G (—£,£). Moreover, since x\ (0) = —w, we have x'(0) + v G {0} x X 2 = T P M. □ 

3 An extension of Theorem 1 and its proof 

In this section we state a version of Theorem 1 extended to non-self-adjoint linear elliptic 
operators of second order and present its proof. 

For that purpose, consider a fixed second order elliptic operator £ 0 under the form 

n n 

£° : = ~ b °i( x ) d i ~ c °( x ) 

i,j =1 i =1 

with coefficients a°-, b ? and c° in (7 r+1 (M n ), r > 0. Note that any second order linear differential 
operator with coefficients in C 1 (M Tl ) can always be placed into the above format. 

Let f^o be a bounded open subset of M” and X and F ,1 (h2 0 ) be as defined in the introduction. 
Given a bounded open subset hi of M n , we denote by A(£ 0 ,f2) C M x ih ( J (O) the set of pairs 
(A q,u), where u is nonzero, satisfying 

f —Cnu = A qu in hi, 

( 2 ) 

u — 0 on dfl 
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Theorem 3. Let f2 0 be a bounded open subset of satisfying the divergence theorem for an 
unit normal vector field u defined on <9f2 0 in almost everywhere and oriented outwards f2 0 . Let 
ip : M —> E 1 (D 0 ) be a C r+1 curve, r > 1, such that </?(0) = la 0 - Denote Tp = <^'(0). Let Ao be 
a multiple eigenvalue of the operator Co and of its adjoint C * 0 with eigenspaces of same 

dimension being generated by, respectively, {wo, u \,..., Uk} and {w 0 , v±,..., Vk} in Hq(Q 0 ) with 
k > 1. Consider the {k + 1) x (k + 1) real matrices A and B with elements, respectively, 


and 


Aij [Ui,Vj) 



a(x) 


dui dvj 
du dv 


<p ■ v dS 


Btj = (ui,Vj ) L 2 := / UiVj dx , 

J S7o 

where a(x ) = j a ij( x ) u i l/ j- Assume that A 0 has algebraic and geometric multiplicities equal 
to k and that the polynomial y(s) = det(A — sB ) has a simple real zero p. Assume also that 
the eigenfunctions Ui and Vi are of C 1 class up to the boundary. Then, there exist C r functions 
t G (—£,£) i ->- A (t) and t G (—£,£) u(t) such that A(0) = Ao, w(0) G (uq, u±, ..., Uk) and 
(A (t),u(t)) G A(—£ 0) £t(^o)) fort G (—£,£). Moreover, we have A'(0) = —p. 

Proof of Theorem 3. Denoted = IBLxHq (Q 0 ) and Y = J8LxH~ 1 (fl 0 ). Let f : E 1 (Q 0 )xX -^-Y 
be the function given by 


f(<p,\,u) = ( - / u 2 dx,T£ o (u ) - A I* Q (u) ) , 
z Jn o 


where Tfi o (u)(v) = B ip( n 0 )((p*u,ip*v), with <p*u(x) = u(<p 1 (a;)) and 


B v {n 0 )(u,v) = LX aij{x)diudjV + bi(x)diuv + c(x)uv dx 

■M^o) itj= i 


2=1 


and WW = J v{ n o) (<p*u, ip*v), with 


J v {n 0 )L, v ) = / uv dx . 

J ip(Oo) 

Note that Tf lo (u) and Iq 0 (u) are linear and continuous, so that / is well-defined. 

Let w G M fc+1 \ {0} be such that A T w = pB T w and consider u* = uuui- Consider the 
function F : M. x X —> Y defined by F(t, A ,u) = f(<p(t), A, it). We now prove that F is in the 
conditions of Theorem [2] at the point (Aq,u*). 



Firstly, F is of C r+1 class by Proposition 2.1 of [8]. Set K = Ker(£ 0 + Ao) C Hq(Q q ) and 
L = Ker(£g + A 0 ) C Hq(Q 0 ). It is clear that 

Im(F> x F(0,A 0 ,i£)) =lx ((«*) © A) CRx H~\n 0 ) 


and 


Ker(Z> x F(0, A 0 , u*)) = {0} x (K D (n*) ± ) CRx //,](O 0 ), 

where denotes the orthogonal in iL _1 (0 0 ) of the subspace M of Hq(Q 0 ). 

Notice that F( 0, Ao, it) = (|, 0) for all it G A' with \\u \\ L 2 = 1. Then, F(0, •) is constant in a 
fc-dirnensional sphere with tangent space at u* equal to {0} x (F fl(it*)- L ) = Ker(F>xF(0, Ao, it*). 
Therefore, the conditions (a) and (b) of Theorem [2] are fulhlled. 

By direct computation, we have 


du 

Ft{ 0, Ao, u) = (0, — (£o + Ao)^=) 

for all u G K. On the other hand, the divergence theorem provides 


du f du dv 

(—(A) + Ao) vp, v) L 2 = / a{x)^^Tp ■ V ds 

(Fp Jdn 0 

for all u G Ker(£ 0 + Ao), v G Ker(£g + Ao). So, writing F t = (F t ! , F 2 ), for any uGf, we derive 


(F t 2 (0, A 0 ,u),n) L 2 = (u,w)^. 


Let u = F 2 (0, Aq, it*) = — (£o + Aq)^=-- By the definition of it*, for each j, we have 


{ll^Vj^jJ2 (U ^ y AjjWj ft ^ , BjjWj , Vj ) , 

meaning that ft — /iit*_Li; for all v E L and u G (u*) © A. 

Then, 


F f (0, A 0 , it*) = (0, it) G Im(F> x F(0, A 0 , it*)), 

so Ft(0,Ao,it*) = D\F(0, Ao, it*)[A, u] for some A G M and u G Hq(Q 0 ), and actually A — fi. 
Thus, the condition (c) of Theorem [2] is satisfied. 
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In order to verify the condition (d) of Theorem [21 we first write the directional derivative 


F^,v){0,X,u) = D x F(0, A, w)[(/q v)} = ((v,u) L 2 , (£ 0 + X)(v) + y,u) 


Thus, for any u G K, 


D(F t - T( M] jj))(0, A 0 , u*)[0, it] = (~{v,u) L 2 , ~(C 0 + A 0 )^= - l^u) 

Notice now that 

coker(iAF(0, A 0 , u*) = R x (#o W)* = ( L n ( U *)-L)* 

v v ' lx ((u*) 0 L x ) 1 ; v x 7 ' 

where the last identification is the natural one. 

Then we need to check that the bilinear function 

(u*)^) x (L n (u*) ± ) ->• M 

G(u,v) = (D(F t 2 - F^)(0, X 0 ,u*)[(0,u)],v) = (u,v) w -n(u,v) 

is non-degenerate. To this end, take a basis of K of the form {uq, u\, , u^}, where u* Q = u* 

and u*-Lu* for i = and a basis of L of the form {uq, v\, . .., vl}, where v*±.u* for 

i — 1,... ,k. Since the algebraic and geometric multiplicities of A 0 are equal, we have u* ^ L 2 -, 
so dim (if n (w*) -1 ) = dim((it*) 0 L 2 -) = k. For i/£l, consider the (k + 1) x (k + 1) matrix 

M{v\, = (ul v*) ¥ - u(ul v*) L 2 = C T ■ (A — uB) ■ D , 

where C and D are the matrices of change of basis. Then 

V 0 N{y) ) 

where N([i) is the k x k matrix of the bilinear form G in the bases ..., ul} and {u*,..., vl}, 
respectively. Then, we have 

det (M(u)) = det(C) det(D) det (A — uB) = (u*, Vq)i?(h — v) det (N(u)) 

and, by the simplicity assumption of ji, we obtain det(lV(/i)) ^ 0. 

So, it follows the first conclusion by Theorem 21 Finally, using the fact that (A'(0), w'(0)) + 
(n,u) G {0} x K, we deduce the desired formula A'(0) = — fi. □ 
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Remark 2. The previous proof can be easily adapted to cover the case of complex zeros ofx(s). 
Precisely, let A 0 be a real eigenvalue of multiplicity k + 1 of the real operator £ 0 and assume 
that the polynomial y(s) has a simple complex zero p. Then, considering the same function f, 
we deduce the existence of a complex branch of eigenvalues A : (—e,e) —$■ C with A(0) = Ao and 
A'(0) = -p. 

4 Applications 

This last section is devoted to some applications to the C r dependence problem of Dirich- 
let eigenvalues with respect to one-parameter C r+1 perturbations of the domain with r > 1. 
Particularly, we consider the question on disks, squares and balls. 

(I) Regularity on disks. Let Do be the unit disk centered at the origin in M 2 and p : R — > 
£d(D 0 ) be a C r+1 one-parameter perturbation such that <p(0) = 1 q 0 . Denote Tp = p'( 0). 

The multiple Dirichlet eigenspaces of the Laplace operator on Do are spanned by functions 
of the form u(p, 9) = f(p) sin (k9) and v(p, 9) = f(p) cos (k9) written in polar coordinates. Here 
fcGN and f(p) = Jk{jp), where Jfc is the Bessel of the first kind of order k and j is some zero 
of J fc . 

Since u and v are L 2 orthogonal, we may assume that / is normalized in such a way that 
B is the identity matrix. The normal derivatives are given by 

du 

— = /'(1) sin (M) 

and 


Let 


dv 

du 


/'(l) cos (kd ). 


5(9) = (p(cos9, sind), (cos 9, sin 9)). 


Let A be the matrix of Theorem |T] Then the discriminant of its characteristic polynomial 
is 


4 b 2 + a 2 


= 4 


- 27 r 


sin (kO) cos(k9)5(9)d9 j 


»2tt 


( cos(k9 ) 2 — sin (k0) 2 )5(9)d9 
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2i r \ ^ / /»27r \ ^ 

sin(2/t6>)5(^)cZ6»J + M cos( 2 fc 0 )<$( 0 )d 0 j = |<f( 2£;)| 2 

which is the 2 /e-Fourier coefficient of S, where a = \ An — A 22 I and b = A± 2 = A 21 . 

Therefore, Theorem [1] provides that all Dirichlet double eigenvalues of the Laplace operator 
on disks are of C r class close to t — 0 whenever 5 ( 2 k) 7 ^ 0 for all k. On the other hand, 
the 2/c-Fourier coefficient of <5 to be nonzero for all k is a generic property with respect to 
perturbations p. 

In order to illustrate this application, we consider the explicit perturbation ft(z) — z + 
where z £ C, (kajJ) G / 1 (M) and a 2 k+ 1 7 ^ 0 for all k. Then, all Dirichlet eigenvalues 
of the Laplace operator on the unit disk are of C°° class close to t — 0. 



(II) Regularity on squares. Let Do = [ 0 , 7 r ] 2 be a square in M 2 and ip : R — > if 1 (Do) be a 
C r+1 one-parameter perturbation such that 99(0) = ln 0 . Denote Tp = ip'(0). 

Multiple eigenspaces of Do with eigenvalue A are spanned by functions of the form 

2 

u a (x i,x 2 ) = - sin (a\Xi) sin (a 2 x 2 ) 

7 T 

with cr 2 + a\ = A and a* G N. 

Let us assume for simplicity that the above Diophantine equation for cr, has only two 
solutions, namely a = (<7i, cr 2 ) and a = (er 2 , cxi), so the eigenspaces are two-dimensional. 
Consider the functions 77, p : [ 0 , 27 t] —> M. 

f ~^ 2 M) if t < 7T f ~<Pi(0,t) if t < 7T 

<y 9 2 (L 7 r) if f > 7 T I p x ( 7 T, f) if t > 7 T 

and the Fourier’s coefficients 


/*27T p2n 

fj(k) = / 77(f) cos (kt)dt , //(&) = / 77(f) cos (kt)dt. 

Jo Jo 

Straightforward computations show that for 7, A G {a, a}, we have 


A. 


■76 


= 1 [ p(t) sin(7 2 f) sin (S 2 t) + (-l) 71+<5l /i( 7r + £) sin(7 2 f) sin(A 2 f)df 

do 


4y 2 


7T 


77(f) sin(7if) sin(A 1 f) + (—1) 72+<S2 77(7 t + f) sin(7 1 f) sin((f 1 f)df 
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and noting that (—l) 71+<Sl — ^_x) 72+<52 , we may change variables in the second term of each 
integrand to obtain 

47-1 Si f 2lT A'yofi? 

A jS = — / //(£) sin( 7 2 t) sm(8 2 t)dt 4-— / ?y(f) sin( 7 if) sin(<5if)df . 

71 J o 71 Jo 

Now using standard trigonometric identities, we have 

A l6 = 2 - - A (72 + $ 2 )) + —^(f?(7i - <*i) - 7(71 + <*i)) • 

7T" 7T Z 

Then the coefficients of A are 

9 a 2 ‘In 2 

A a a = -4(A(0) - £(2a 2 )) + -^W(O) - rj(2a 1 ))), 

7T Z 7T Z 

2rr 2 2rr 2 

= —^(A(0) - £(2^)) + -^W(O) - 7}(2a 2 ))), 

7T Z 7T Z 

Arif = A aW = ~ Gl ° 2 (£((7i - cr 2 ) - £(<71 + cr 2 ) + 7)0l - cr 2 ) - 7(07 + cr 2 )) • 

7T Z 

Finally, A has simple eigenvalues if and only if either A aa 7 ^ A-^ or A a ^ ^ 0. Or equivalently, 
if and only if one of the following conditions occur: 

(v + £)(oi + ^2) 7^ (f) + £)Oi - *2) , 

0 ? (£ - 7))(2cr 2 ) - of(£ - 1 ?)(2c r i) ^ (<r? - (T 2 )(£ - fj)( 0) . 

Thus, since these two situations are generic with respect to perturbations p, Theorem [T| pro¬ 
vides that all Dirichlet eigenvalues of the Laplace operator on squares are of C r class close to 
t = 0 for almost every perturbation p. 

(Ill) Regularity on balls in M n . Let Oo be the unit ball centered at the origin in M n 
and p : R —* £ ,1 (f2 0 ) be a C r+1 one-parameter perturbation such that <^(0) = ln 0 . Denote 

p = p'(0). 

The first multiple eigenspace in the unit ball of M n is spanned by functions of the form 
Ui(x) = where p = ||x||, i = 1,.. .,n. The derivatives in the normal direction coincide 

in S n_1 with aq,..., x n , except for a constant. 
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Denote 


8{0) = (&{0),v{Q)), 

where v is the unit normal vector field of § n_1 oriented outwards. Its Fourier coefficients of 
degree 2 are 


Fi= / (et-9i)5(0)dS,i = 2,...,n, 

J S "- 1 


Cij = / 9idj5(9)dS , = 

Let a = Jgn-i &{5{6)dS, we compute 

( 


A — aid = 


V 


0 

C 21 

C 31 . 

.. C nl 

C'2l 

f 2 

C 32 ■ 

.. C n2 

•• p 

C 32 

F 3 . 

.. Cn 2 

c n 1 

C n 2 

Cn 3 .. 

• • F n 


\ 




Since Xa( s ) = XA-aid(s ~ a), the simplicity of the eigenvalues of A can be translated from the 
second order Fourier coefficients of 5. In particular, the characteristic polynomial xa of the 
matrix A is constructed from sums and products of these latter ones. Thus, the discriminant 
of Xa is nonzero for almost every perturbation <p, so that all n eigenvalues are of C r class close 
to t — 0. 

The problem of determining the simplicity of the eigenvalues of A in general seems to be a 
complicated one for balls in M n with n > 3 and we don’t hope to find a simple solution in these 
cases. 


(IV) Two disjoint domains. Let us apply Theorem Q] into a quite situation. 

Consider C r+1 one-parameter perturbations ip : M —» I? 1 (Do) and i/j : R —> I? 1 (Vo)- Denote 
U t = Pt{U 0 ) and V t = iptiYo) with principal eigenvalues and u t , respectively, where U 0 and 
Vo are disjoint bounded open subsets of M n but equal after a translation. 

Let W — Uq U Vo- Since p 0 = the eigenspace is spanned by functions u t and v t with 
support at U t and V tl respectively. Assuming that ||wo||l 2 = II^oIIl 2 = 1, the inner product 
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matrix B is the identity and A is diagonal. Moreover, An = —/i'( 0) and A 2 2 = — i/(0) by the 
Hadamard formula, see for example [Sj. 

In this case we evoke Theorem Q] which states that (it and u t are C r functions on t, provided 
that kin ~f~ A 2 2 - However, observe that, unlike u t and v t , not every eigenfunction of W t = U t CV t 
is smooth on t. For example, if wo = ug + vq the local uniqueness of the principal eigenvalues 
of U 0 and V 0 guarantees that the unique possible continuation is w t = u t + v t which is not an 
eigenfunction unless gt t = v t for every t small enough. Part of the proof of Theorem [I] is to find 
the eigenfunction having a continuation. 
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